Astronomy Sz Astrophysics manuscript no. CoRoT2'IRAC 


November 26, 2009 


(DOI: will be inserted by hand later) 





The thermal emission of the young and massive planet 
CoRoT-2b at 4.5 and 8 /xm* ** 

M. Gillon 1 ' 2 , A. A. Lanotte 1 , T. Barman 3 , N. Miller 4 , B.-O. Demory 2 , M. Deleuil 5 , J. Montalban 1 , 
F. Bouchy 6 , A. Collier Cameron 7 , H. J. Deeg 8 , J. J. Fortney 4 , M. Fridlund 9 , J. Harrington 10 , P. Magain 1 , 
C. Moutou 5 , D. Queloz 2 , H. Rauer 11 ' 12 , D. Rouan 13 , J. Schneider 14 

1 Institut d'Astrophysique et de Geophysique, Universite de Liege, Allee du 6 Aout 17, Bat. B5C, 4000 Liege, 
Belgium 

2 Observatoire de Geneve, Universite de Geneve, 51 Chemin des Maillettes, 1290 Sauverny, Switzerland 

3 Lowell Observatory, 1400 West Mars Hill Road, Flagstaff, AZ 86001, USA 

4 Department of Astronomy and Astrophysics, University of California, Santa Cruz, USA 

5 LAM, UMR 6110 CNRS, 38 rue Frederic Joliot-Curie, 13388 Marseille, France 

6 Observatoire de Haute Provence, USR 2207 CNRS, OAMP, F-04870 St-Michel l'Observatoire, France 

7 School of Physics and Astronomy, University of St. Andrews, North Haugh, Fife, KY16 9SS, UK 

8 Instituto de Astroh'sica de Canarias, C. Via Lactea S/N, E-38200 La Laguna, Spain 

9 Research and Scientific Support Department, European Space Agency, ESTEC, 220 Noordwijk, The Netherlands 

10 Planetary Sciences Group, Department of Physics, University of Central Florida, Orlando, FL 32816, USA 

11 Institut fuer Planetenforschung, DLR, Rutherford str. 2,12489 Berlin, Germany 

12 Zentrum fuer Astronomie und Astrophysik, Hardenbergstr. 36,10623 Berlin, Germany 

13 LESIA, UMR 8109 CNRS, Observatoire de Paris, UVSQ, Universite Paris-Diderot, 5 Place J. Janssen, 92195 
Meudon, France 

14 LUTH, UMR 8102 CNRS, Observatoire de Paris-Meudon, 5 Place J. Janssen, 92195 Meudon, France 
Received date / accepted date 

Abstract. We report measurements of the thermal emission of the young and massive planet CoRoT-2b at 4.5 and 
8 /jm with the Spitzer Infrared Array Camera (IRAC). Our measured occultation depths are 0.510±0.042 % and 
0.41 ± 0.11 % at 4.5 and 8 fim, respectively. In addition to the CoRoT optical measurements, these planet/star 
flux ratios indicate a poor heat distribution to the night side of the planet and are in better agreement with an 
atmosphere free of temperature inversion layer. Still, the presence of such an inversion is not definitely ruled out 
by the observations and a larger wavelength coverage is required to remove the current ambiguity. Our global 
analysis of CoRoT, Spitzer and ground-based data confirms the large mass and size of the planet with slightly 
revised values (M p = 3.47 ±0.22 Mj, R p = 1.466 ±0.044 Rj). We find a small but significant offset in the timing 
of the occultation when compared to a purely circular orbital solution, leading to ecosaj = —0.00291 ± 0.00063 
where e is the orbital eccentricity and u is the argument of periastron. Constraining the age of the system to 
be at most of a few hundreds of Myr and assuming that the non-zero orbital eccentricity is not due to a third 
undetected body, we model the coupled orbital-tidal evolution of the system with various tidal Q values, core sizes 
and initial orbital parameters. For Q' a — 10 5 — 10 6 , our modelling is able to explain the large radius of CoRoT-2b 
if Q'p < 10 J ' 5 through a transient tidal circularization and corresponding planet tidal heating event. Under this 
model, the planet will reach its Roche limit within 20 Myr at most. 
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* Based on data collected with the VLT/FORS2 instrument 
at ESO Paranal Observatory, Chile (programs 081.C-0413(B)). 

** The photometric time-series used in this work arc 
only available in electronic form at the CDS via anony- 
mous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
|http://cdsweb.u-str asbg.fr/cgi-bin/qcat ?J/A±A/| 



1. Introduction 

Transiting planets are key objects for our understanding 
of the atmospheric properties of exoplanets. Indeed, their 
special geometrical configuration gives us the opportunity 
not only to deduce their density but also to study di- 
rectly their atmospheres without the challenging need to 
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spatially resolve their light from that of their host star. 
In particular, their emergent flux can be directly mea- 
sured during their occultation (secondary eclipse) when 
they are hidden by their host star, as was demonstrated 
by Charbonneau et al. (2005) and Deming et al. (2005). 
Since 2005, many exoplanet occultation measurements 
have been gathered, the bulk of them by the Spitzer 
Space Telescope (see, e.g., Deming 2009), the few oth- 
ers being due to the Hubble Space Telescope (Swain et 
al. 2009), CoRoT (Alonso et al. 2009b, 2009c; Snellen 
et al. 2009a), Kepler (Borucki et al. 2009) and ground- 
based telescopes (Sing & Lopez-Morales 2009, de Mooij 
& Snellen 2009, Gillon et al. 2009b). Combining the pho- 
tometric measurements at different wavelengths allows us 
to map the Spectral Energy Distribution (SED) of the 
planet and to constrain its chemical composition, its ther- 
mal distribution efficiency and the presence of a possible 
stratospheric thermal inversion (see, e.g., Charbonneau et 
al. 2008, Knutson et al. 2008). Such inversions have been 
detected for the highly irradiated planets HD 209458b 
(Burrows et al. 2007b, Knutson et al. 2008), TrES-2b 
(O'Donovan et al. 2009), TrES-4b (Knutson et al. 2009), 
XO-lb (Machalek et al 2008) and XO-2b (Machalek et al. 
2009). These results are in rather good agreement with 
the theoretical division of hot Jupiters into two classes 
based on their level of irradiation (Hubeny et al. 2003, 
Burrows et al. 2007b, Harrington et al. 2007, Burrows et 
al. 2008, Fortney et al. 2008). Under this division, the 
planets warmer than required for condensation of high- 
opacity gaseous molecules like TiO/VO, tholins or poly- 
acetylenes should show a stratospheric temperature inver- 
sion due to the absorption in their upper-atmosphere of 
a significant fraction of the large incident flux by these 
compounds. The less irradiated planets would lack these 
gazeous compounds and the resulting temperature inver- 
sion. Still, this simple division was recently challenged by 
the absence of thermal inversion reported for the strongly 
irradiated planet TrES-3b by Fressin et al. (2009). This 
result indicates that, in addition to the irradiation am- 
plitude, other effects like chemical composition, surface 
gravity and the stellar spectrum have probably an impact 
on the temperature profile of hot Jupiters. 

With an irradiation ~ 1.3 x 10 9 erg s _1 cm~ 2 , the 
planet CoRoT-2b could be expected to show such a tem- 
perature inversion, according to the theoretical division 
mentionned above. This planet is the second one discov- 
ered by the CoRoT transit survey mission (Alonso et al. 
2008, hereafter A08). Spectral analysis and evolution mod- 
elling of the host star leads to a solar-type dwarf with a 
mass M* = 0.97 ± 0.06 M Q and an effective temperature 
T cff = 5625 ± 120 K (Bouchy et al. 2008, hereafter B08). 
A08 derived for the planet a radius of 1.465 ± 0.029 Rj 
and a mass of 3.31 ± 0.16 Mj, leading to a density of 
1.31 ± 0.04 g cm~ 3 , very close to the value for Jupiter. 
This density is surprising because the radius of massive 
planets is expected to approach Jupiter's asymptotically. 
In this context, it is worth noticing the probably young 
age of the system. Indeed, the presence in the stellar spec- 



trum of the Li I absorption line and the strong emission 
in the Ca II H and K line cores (B08) suggest that the 
star is still close to the Zero- Age Main-Sequence (ZAMS) 
and is thus younger than 0.5 Gyr (B08), in full agreement 
with the short rotational period of ~ 4.5 days deduced 
from CoRoT photometry (A08). Still, the youngness of 
CoRoT-2b is not enough to prevent it from falling into 
the sub-group of planets with a radius larger than pre- 
dicted by basic models of irradiated planets (Burrows et 
al. 2007a, Fortney et al. 2007). Most of these planets show 
an orbital eccentricity compatible with zero. Nevertheless, 
these planets could still have undergone during their evo- 
lution a tidal heating large enough to explain their low 
density (Jackson et al. 2008b; Ibgui & Burrows 2009), and 
it is thus important to measure very precisely their present 
eccentricity to constrain their tidal and thermal history. 
The precise measurement of a planet's occultation pro- 
vides strong constraints on the orbital eccentricity, espe- 
cially on the parameter ecosu->, where e is the eccentricity 
and u> the argument of periastron (see e.g. Charbonneau 
et al. 2005, Knutson et al. 2009). In the case of CoRoT-2b, 
the dynamical interest of such occultation observations is 
reinforced by the large jitter noise of the young host star 
(B08), which makes the precise determination of a tiny 
eccentricity very challenging for the radial velocity (RV) 
method alone. 

With the goals of better characterizing the atmo- 
spheric properties of CoRoT-2b (SED, inversion) and im- 
prove our understanding of its low density (tidal heating) , 
we observed the occultation of this planet at 4.5 and 8 
/im with Spitzer /IRAC (DDT program 486). A partial 
transit was also observed with VLT/FORS2 to bring one 
more constraint on the orbital parameters. We report here 
the results of the analysis of these new data. Section 2 
presents our IRAC and VLT observations and their re- 
duction. We analyzed this new photometry in combina- 
tion with CoRoT transit photometry and published RVs. 
This combined analysis is presented in Sec. 3. We present 
and discuss our results in Sec. 4 and give our conclusions 
in Sec. 5. 

2. New photometric observations 

2.1. IRAC occultation photometry 

CoRoT-2 (2MASS 19270649+0123013, K s = 10.31) was 
observed by Spitzer (Werner et al. 2004) during an oc- 
cultation of its planet on November I s * 2008 from 03h50 
to 08h50 UT. The observations were performed with the 
Infrared Array Camera (IRAC) (Fazio et al. 2004) in full 
array mode (256 x 256 pixels, 1.2 arcsec/pixel) simultane- 
ously at 4.5 and 8 /im. The telescope was not repointed 
during the course of the observations to minimize the 
motion of the stars on the array. We carefully selected 
the pointing in order (1) to avoid the bright star 2MASS 
19270954+0123280 {K s = 7.55) that would have saturated 
the detector for any exposure time while ensuring that it 
will not fall into one of several regions outside the FOV 
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that arc known to result in significant scattered light on 
the detectors, and (2) to avoid areas of the array with 
known bad pixels or significant gradients in the flat field 
as well as areas known to be affected by scattered starlight. 
We ensured also that no bright star would have been lo- 
cated into stray light avoidance zonefl An effective in- 
tegration time of 10.4s was used during the whole run, 
resulting in 1385 images for each channel. For our analy- 
sis, we used the images calibrated by the standard Spitzer 
pipeline (version S18.0) and delivered to the community as 
Basic Calibrated Data (BCD). We converted fluxes from 
the Spitzer units of specific intensity (MJy/sr) to pho- 
ton counts, and aperture photometry was obtained for 
CoRoT-2 in each image using the IRAF/DAOPHOT0 soft- 
ware (Stetson, 1987). In both channels, the Point-Spread 
Function (PSF) of the target is slightly blended with the 
one of the fainter (K a = 12.03) redder (J - K s = 0.84 
vs 0.47 for CoRoT-2) star 2MASS 19270636+0122577 lo- 
cated at ~ 4" (see Fig. 1). A small aperture radius was 
used for both channels (4.5 fim: 4 pixels, 8 /zm: 3.5 pix- 
els). The aperture was centered in each image by fitting 
a Gaussian profile on CoRoT-2. A mean sky background 
was measured in an annulus extending from 8 to 16 pix- 
els from the center of the aperture, and subtracted from 
the measured flux for each image. Each measurement was 
compared to the median of the ten adjacent images and 
rejected as an outlier if the difference was larger than four 
times its theoretical error bar. Twenty-four points (3.5 %) 
were rejected at 4.5 /xm and 37 points (2.7 %) at 8 fim. 
Figure 2 shows the resulting time-series for both channels. 
Despite its small size, the photometric aperture does not 
only contain counts due to CoRoT-2 but also to the nearby 
fainter star, leading to a dilution of the eclipse. To esti- 
mate this dilution and correct the measured eclipse depths 
for it, we used the following procedure. For both channels, 
we partially deconvolved the images taken after the occul- 
tation, using the deconvolution program DECPH0T (Gillon 
et al. 2006, 2007a, Magain et al. 2007). We used the over- 
sampled high-SNR PSF available on Spitzer's web sited to 
deduce the partial PSF needed for this deconvolution. The 
deconvolved images (see Fig. 1) are oversampled by a fac- 
tor of 2 and their PSF is a Gaussian with a Full- Width at 
Half Maximum (FWHM) of 2 pixels, corresponding thus 
to a FWHM of 1 pixel for the original sampling. This has 
to be compared to a FWHM ~ 1.5 pixel for the original 
images. We performed aperture photometry on the PSF 
model to measure the fraction a of the flux of CoRoT-2 
within an aperture of 8 (4.5 fim) and 7 (8 fim) pixels. At 
this stage, we compared the total flux F to t of CoRoT-2 
obtained by DECPH0T for each image to the flux obtained 
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For details, see the IRAC Data Handbook available at 
|http:/ /ssc. spitzer.caltech.edu/irac] 

1 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under cooper- 
ative agreement with the National Science Foundation. 

3 http:/ /ssc. spitzer. caltech.edu/irac/psf.html 





Fig. 1. Top: Zoom on CoRoT-2 and the nearby fainter 
star within an IRAC image taken at 4.5 /jm (left) and 8 
/im (right). Bottom: Same, but after deconvolution. 



with aperture photometry on the model images (i.e. the 
obtained higher-resolution images convolved by the partial 
PSF model). This last measurement should be the sum of 
a x F tot and the contaminating flux due to the nearby star. 
Subtracting a x F tot to this quantity and dividing by the 
same a x F tot finally gave us an estimation of the aperture 
contamination due to the nearby star. Considering all the 
images taken after the occultation, we obtained a dilution 
of 16.4 ± 0.4 % and 14.3 ± 0.7 %, respectively at 4.5 and 
8 fim. 

2.2. VLT/FORS2 transit photometry 

A partial transit of CoRoT-2b was observed on September 
9 th 2008 with the FORS2 camera (Appenzeller et al. 1998) 
installed at the VLT/UT1 (Antu). The FORS2 camera has 
a mosaic of two 2k x 4k MIT CCDs and is optimized for 
observations in the red with a very low level of fringes. It 
was used several times in the past to obtain high preci- 
sion transit photometry (e.g. Gillon et al. 2007b, 2009a, 
2009b; Pont et al. 2007). The high resolution mode and 
lxl binning were used to optimize the spatial sampling, 
resulting in a 4.6' x 4.6' field of view with a pixel scale 
of 0.063" /pixel. Airmass decreased from 1.18 to 1.11 then 
increased to 1.35 during the run which lasted from 23h40 
to 3hl2 UT. The quality of the night was photometric. 

We acquired 448 images in the z-GUNN+78 filter 
(X eff = 910 nm, FWHM = 130.5 nm) with an expo- 
sure time ranging from 0.6 to 3 s. After a standard pre- 
reduction, the stellar fluxes were extracted for all the im- 
ages with the IRAF/DA0PH0T aperture photometry soft- 
ware. Fifty images were revealed to be saturated and were 
discarded from the analysis. Several sets of reduction pa- 
rameters were tested, and we kept the one giving the most 
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Fig. 2. IRAC occultation photometry obtained at 4.5 /xm 
(top) and 8 /im (bottom). For both time-series, the best- 
fitting occultation+systematics model is superimposed (in 
red). 

precise photometry for the stars of similar brightness than 
CoRoT-2. After a careful selection of reference stars, dif- 
ferential photometry was obtained. The resulting transit 
light curve is shown in Fig. 3. After subtraction of the 
best-fit model (see next section), the obtained residuals 
show a standard deviation of ~ 1.9 x 10 -3 , close to the 
mean theoretical noise (~ 1.7 x 10~ 3 ). 

3. Data analysis 

3.1. Used data and model 

We performed a global determination of the system pa- 
rameters based on our new photometry in addition to the 
following data: 

— The phase-folded CoRoT transit photometry pre- 
sented in A08. The 160 measurements of this transit 
light curve were obtained after folding the 78 tran- 
sits observed by CoRoT using the precise ephemeris 
deduced in A08 and after binning the resulting light 
curve with a bin size of ~ 2.5 min. For our analy- 
sis, we projected this phase-folded photometry onto 
the central transit timing presented in A08. To take 
into account the uncertainty on the time of mini- 
mum light and on the orbital period presented in 
A08 (T = 2454237.53562 ± 0.00014 HJD and P = 
1.7429964 ± 0.0000017 days), new values were ran- 
domly drawn from the corresponding normal distri- 
butions at the beginning of each chain of the MCMC 



Fig. 3. Top: VLT/FORS2 z-band transit photometry with 
the best-fitting transit+systematics model superimposed 
(in red). Middle and bottom: Residuals of the fit unbinncd 
and binned per 20 minutes. The larger scatter during the 
transit is probably due to the inhomogeneity of the stellar 
surface (spots). 

analysis (see below) before projecting the phase-folded 
light curve. 

— The radial velocity (RV) measurements published in 
A08 and B08 and obtained by the HARPS and 
SOPHIE spectrographs. These RVs encompass two 
transits. These spectroscopic transit observations were 
obtained to measure the sky-projected angle /3 between 
the planetary orbital axis and the stellar rotation axis 
via the observation of the Rossiter-McLaughlin effect 
(RM; Queloz et al. 2000). We included these spectro- 
scopic transit observations in our analysis to benefit 
from as many constraints as possible on the orbital 
and eclipses parameters. 

These data were used as input into an adaptative 
Markov Chain Monte Carlo (MCMC; see e.g. Tegmark 
2004, Gregory 2005, Ford 2006) algorithm. MCMC is a 
Bayesian inference method based on stochastic simula- 
tions that samples the posterior probability distribution of 
adjusted parameters for a given model. Our MCMC imple- 
mentation uses the Metropolis-Hasting algorithm (see e.g. 
Carlin & Louis 2008) to perform this sampling. Our model 
is based on a star and a transiting planet on a Keplerian 
orbit about their center of mass. More specifically, we used 
a classical Keplerian model for the RVs obtained outside 
the transit in addition to a Rossiter-McLaughlin effect 
model (Gimenez 2006) for the RVs obtained during tran- 
sit. To model the eclipse photometry, we used the photo- 
metric eclipse model of Mandel & Agol (2002) multiplied 
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by a systematic effect model different for each time-series 
(see Sec. 3.3). 

3.2. Limb-darkening 

For both photometric transits, a quadratic limb-darkening 
law was assumed. For the FORS2 light curve, the 
quadratic coefficients u\ and u 2 were kept fixed to 0.23 and 
0.32, the values deduced from Claret's tables (2000; 2004) 
for T eff = 5625 K, log g = 4.3 and [Fe/H] = 0.0 (B08). 
Considering the excellent quality of the CoRoT transit 
photometry, we allowed for it the quadratic coefficients 
U\ and u 2 to float in our MCMC analysis, using as jump 
parameter^ not these coefficients themselves but the com- 
binations ci = 2 x u\ + u 2 and c 2 — u\ — 2 x u 2 to minimize 
the correlation of the obtained uncertainties (Holman et 
al. 2006). To obtain a limb-darkening solution consistent 
with theory, we decided to use a Bayesian penalty on c\ 
and c 2 based on theoretical values and errors for u\ and 
u 2 . The broad CoRoT band-pass does not correspond to 
any photometric filter, but its maximum of transmission is 
close to the V and R bands (Deleuil et al. 2008). We used 
the method described in Gillon et al. (2009b) to deduce 
from Claret's tables (2000; 2004) the theoretical values for 
ui and u 2 and their errors <j Ul and <j U2 for the V and R 
filters and the spectroscopic parameters of CoRoT-2b re- 
ported in B08. For each coefficient, we took as initial value 
the mean of the values obtained for both filters. For the 
errors, we took the mean of the errors deduced for both 
filters and added it quadratically to the difference between 
both filters to take into account our ignorance of the ef- 
fective wavelength of the photometry. We obtained this 
way iti = 0.413 ± 0.108 and u 2 = 0.293 ± 0.038 for our 
initial limb-darkening coefficients. Finally, the following 
Bayesian penalty was added to our merit function: 

where c! i is the initial value deduced for the coefficient Ci 
and <7 C ' is its error computed from a Ul and a U2 . 

For the spectroscopic transits, a quadratic limb- 
darkening law was also assumed. The values u\ = 0.465 
and u 2 — 0.276 were deduced from Claret's tables for 
the stellar parameters presented in B08 and for the V- 
filter, corresponding to the maximum of transmission of 
the HARPS and SOPHIE instruments. These values were 
kept fixed in the MCMC. 

3.3. Modelled photometric systematic effects 

For each light curve, the eclipse model was multiplied by 
a trend model to take into account known low-frequency 
noise sources (instrumental and stellar). 

At 4.5 /im (InSb detector), the measured IRAC fluxes 
show a strong correlation with the position of the tar- 

4 Jump parameters are the model parameters that are ran- 
domly perturbed at each step of the MCMC. 
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get star on the array. This effect is due to the inhomo- 
geneous intra-pixel sensitivity of the detector and is now 
well-documented (see, e.g., Knutson et al. 2008 and refer- 
ences therein). Following Charbonneau et al. (2008), we 
modelled this effect with a quadratic function of the sub- 
pixel position of the PSF center: 

A(dx, dy) = ai + a 2 dx + a 3 dx 2 + a^dy + a^dy 2 + a^dxdy ,(2) 

where dx and dy are the distance of the PSF center to the 
center of the pixel. Notice that we followed here Desert 
et al. (2009) and added the cross-term ae to the function 
A(dx,dy). We measured the PSF center for CoRoT-2 in 
each image by fitting a Gaussian profile. Its x position 
ranged from 192.89 to 193.05 during the run, while its y 
position ranged from 240.12 to 240.39. 

At 8 /im (SiAs dectector), the intra-pixel sensitivity 
homogeneity is good, but another systematic affects the 
photometry. This effect is known as the 'ramp' because it 
causes the gain to increase asymptotically over time for 
every pixel, with an amplitude depending on their illumi- 
nation history (see e.g. Knutson et al. 2008 and references 
therein). Following Charbonneau et al. (2008) again, we 
modelled this ramp as a quadratic function of ln(dt): 

B{dt) = bi + b 2 ln{dt) + b 3 (ln(dt)) 2 , (3) 

where dt is the elapsed time since 15 min before the start 
of the run. 

From the CoRoT photometry, the star CoRoT-2 is 
known to be variable at the 2-3% level on a timescale of 
~ 4.5 days, corresponding to its rotational period (A08; 
Lanza et al. 2009). For the VLT and IRAC time-series, 
we modeled this low-frequency modulation by a time- 
dependent quadratic polynomial: 

C(dt) = cx + c 2 dt + c 3 dt 2 , (4) 

where dt is the elapsed time since 15 min before the start 
of the run. As the photometric modulation due to rotating 
spots is a wavelength-dependent effect, independent coef- 
ficients were fitted for the two IRAC time-series despite 
their covering of the same occultation. 

At the end, the VLT/FORS2 trend model thus has 
three coefficients, the IRAC 4.5 fim 3 + 6-1 = 8 coeffi- 
cients and the IRAC 8 /*m 3 + 3-1 = 5 coefficients. All 
the used trend models are linear in their coefficients, so 
instead of considering these coefficients as jump parame- 
ters in the MCMC, we choose to determine them by linear 
least squares minimization at each step of the MCMC af- 
ter division of the data by the eclipse model generated 
from the latest set of jump parameters (see Sec. 3.6). We 
used for this purpose the SVD method (Press et al. 1992), 
which has been found to be very robust. 

The CoRoT transit photometry is already corrected 
for known systematics and the stellar rotational variabil- 
ity. Nevertheless, we preferred to do not assume it to be 
perfectly normalized and consider a flux normalisation fac- 
tor d norm that was also determined via SVD at each step 
of the MCMC. 
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3.4. Photometric correlated noise 

Taking into account the correlation of the noise is impor- 
tant to obtain reliable error bars on the fitted parameters 
(Pont et al. 2006). For this purpose, we followed a proce- 
dure similar to the one described by Winn et al. (2008). 
For each light curve, the standard deviation of the resid- 
uals of the first chain was determined for the best-fitting 
solution, without binning and with several time bins rang- 
ing from 10 to 30 minutes. For each binning, the following 
factor /3 re( j was determined 

where N is the mean number of points in each bin, M is 
the number of bins, and o\ and &n are respectively the 
standard deviation of the unbinned and binned residu- 
als. The largest value obtained with the different binnings 
was used to multiply the error bars of the measurements. 
We obtained /3 rcd = 2.2 for the VLT/FORS2 light curve, 
1.25 for the IRAC/4.5/im curve, 1.14 for the IRAC/S/zm 
curve and 1.25 for the CoRoT photometry Thus, all the 
light curves show a significant level of correlated noise. 
The large /3 re( j obtained for the FORS2 light curve can be 
attributed to the presence of spots on the stellar surface 
and the resulting increase of the correlation of the noise 
during the transit (see Fig. 3). 

3.5. Systemic RVs and jitter noise 

For each RV time-series, the systemic velocity was deter- 
mined at each step of the MCMC from the residuals via 
SVD. Our code is able to account for more linear terms, 
i.e. for trends in the RV time-series, but it was not needed 
here. For the RV data taken outside transit, we assumed 
a different systemic velocity for the SOPHIE and HARPS 
data to account for a possible difference of zero-point 
calibration between both instruments. Following B08, we 
added quadratically a jitter noise of 56 m s _1 to the errors 
to account for the stellar activity. For the spectroscopic 
transit data, we considered the same RV offset during the 
whole transit, so we did not add any jitter noise but only 
considered a different systemic velocity for both spectro- 
graphs. Our analysis of the residuals of the first MCMC 
chain showed us that the jitter noise of 56 m s _1 assumed 
for the data taken outside transit was leading to a residual 
rms in good agreement with the mean error of the mea- 
surements. Still, we had to add an extra-noise of 13 m s _1 
for the data taken during transit to obtain a similar agree- 
ment. The need for this extra-noise could be explained by 
the inhomogenous surface of the spotted star and/or by 
the systematic errors brought by the measurement of the 
RV via the fit of a Gaussian profile on the non symmetric 
cross correlation function of the spectrum (see Winn et al. 
2005, Triaud et al. 2009). 

One could wonder why the correlation of the noise is 
not treated in a similar way for the spectroscopic and 
photometric eclipse time-series. The answer is that the 



time sampling of the RV time-series is much poorer than 
the one of the light curves and is similar to the timescale 
of ingress/egress. We can thus not estimate precisely the 
level of correlated noise at this frequency via the method 
described in Sec. 3.4. Still, the fact that the time sampling 
and the correlation timescale of the noise that we want to 
model are similar makes the addition of the quadratic dif- 
ference between the residual rms and the mean RV error 
a proper method to take into account this 'red' noise. 

3.6. Jump parameters, priors and merit function 

The jump parameters in our MCMC simulation were: 
the planet/star area ratio (R p /R s ) 2 , the transit width 
(from first to last contact) W, the impact parameter 
b' = acosz/i?,, the two Lagrangian parameters ecosw 
and esinw where e is the orbital eccentricity and u is the 
argument of periastron, and the Ki parameter character- 
izing the amplitude of the orbital RV signal (see Gillon et 
al. 2009b). We assumed a uniform prior distribution for 
all these jump parameters. 

The products V sin/ cos (3 and V sin/ sin /? were also 
jump parameters in our MCMC, where Vsinl is the pro- 
jected stellar rotational velocity and ft is the spin-orbit 
angle (see Gimenez 2006). As we have an independent 
determination of Vsinl from spectroscopy (10.5 ± 0.4 
km s _1 , B08), we added the following Bayesian penatly 
to our merit function: 

RP (VsinJ - VsmI B0S ) 2 

°V sin I B os 

where VainlBos is 10.5 km s _1 and cry s in/ BO s i s 0-4 
km s . 

A totally independent determination of the orbital pe- 
riod P and time of minimum light T was impossible be- 
cause we folded the CoRoT transit photometry with the 
ephemeris presented in A08. This is why we let these pa- 
rameters vary under the control of the following Bayesian 
penalty: 

Rp (P-Paos) 2 , (T-T Am f 

ephemeris 2 '2 1 v / 

° F 'aos ^Taos 

where Pao8 and Taos are the best-fitting values presented 
in A08 and ap A0S and <TT Aoa are their errors. In other 
words, we used a normal prior distribution for these two 
jump parameters based on the CoRoT results reported in 
A08. 

The merit function used in our analysis was the sum of 
the x 2 f° r each time-series and of the Bayesian penalties 
presented in Eq. 1, 6 and 7. 

3. 7. Structure of the analysis 

Our analysis was similar to the one presented by Gillon et 
al. (2009b), consisting in four successive steps: 

1. First, we performed a single MCMC chain aiming to 
assess the level of correlated noise in the photometry 
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and of jitter noise in the RVs and to update the mea- 
surement error bars accordingly. This chain was com- 
posed of 10 5 steps, the first 20 % of each chain being 
considered as its burn-in phase and discarded. 
Five new MCMC chains (10 5 steps each) were per- 
formed using the updated measurement error bars. 
The good convergence and mixing of these five chains 
was checked succesfully using the Gelman and Rubin 
(1992) statistic. The inferred value and error bars for 
each parameter were obtained from the marginalized 
posterior distribution. The goal of this second step was 
to provide us with an improved estimation of the stel- 
lar density p* = 1.3lj^§| Pq. 

The deduced stellar density and the spectroscopic pa- 
rameters presented in B08 were then used to determine 
the stellar mass M* and age r* via a comparison with 
the stellar evolution models computed with the CLES 
code (Scuflaire et al. 2008). We obtained a stellar mass 
M* = 0.96±0.08M© and a stellar age r* = 2.7±|;f Gyr. 
A new run of 20 MCMC chains was then performed. 
This step was identical to the second one, with the 
exception that at each step of the chains, the physical 
parameters M p , R p and i?„ were computed from the 
relevant jump parameters and the stellar mass. For this 
latter, a value was randomly drawn at each step from 
the normal distribution 7V(0.96, 0.08 2 )M Q derived in 
the previous step. 



4. Results and discussion 

Table 1 shows the median values and 68.3% probability 
interval for the jump and physical parameters given by 
our MCMC simulation, and compares them to the val- 
ues presented in A08/B08. Notice that the planet/star 
flux ratios reported in Table 1 are the deduced occulta- 
tion depths corrected for the signal dilution due to the 
nearby star (see Sec. 2.1). Figure 4 shows the IRAC pho- 
tometry corrected for the systematic and binned per five 
minutes, with the best-fitting eclipse model superimposed. 
The best-fitting models for the CoRoT and spectroscopic 
transits are presented in Fig. 5. 

4.1. CoRoT-2b: a young, bloated and massive planet 
in a slightly eccentric and well-aligned orbit 

As shown in Table 1, our results for CoRoT-2b agree well 
with the results reported in A08 and B08. Our error bars 
are in average larger than the ones reported in these pre- 
vious works. We consider our error bars as more reliable 
for the following reasons: (1) we did not assume a circular 
orbit, (2) we took into account the ephemeris errors for 
the folded CoRoT photometry, (3) we did not assume a 
perfect normalization for the CoRoT photometry, and (4) 
we took into account the correlated noise present in the 
light curves. To assess the impact of the VLT photometry 
on the final solution, we performed another MCMC inte- 
gration without it, letting P and To free but under the 
control of a Bayesian prior based on A08 ephemeris. The 




4771.75 4771.8 
HJD - 2450000 



Fig. 4. IRAC 4.5 ^m (top) and 8 pm (bottom) occultation 
photometry binned per five minutes and corrected for the 
systematics with the best-fitting occultation models su- 
perimposed. The bottom dataset is shifted for clarity. The 
dilution of the occultation due to the nearby fainter star 
is not corrected here. 




E 




dT (days 



Fig. 5. Top: CoRoT transit photometry with the 
best-fitting transit model superimposed. Bottom: 
HARPS/SOPHIE transit RVs with the best-fitting 
RM model superimposed. 



resulting parameters and their error bars were all in good 
agreement with the values shown in Table 1, the main ef- 
fect of the VLT light curve being to improve by a factor 
~ 2 the precision on the orbital period because it was ob- 
tained nearly one year after the end of the last CoRoT 
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Parameter 



This study 



A08/B08 



Unit 



Jump parameters 



Planet/star area ratio (R p /R s ) 2 

b' = a cos i/R* 

Transit width W 

e cos lo 

esinw 

V sin / cos (5 

V sin / sin f3 
4.5 fim dF 2 
8 /xm dF?. 
RV K 2 



0.02750 ± 0.00012 

r, 99 O + 0.018 
U.^^O_g Q2Q 

0.09446i° :°~i; 
-0.0029li°« 
0.01391°;°°™ 
10.79l°;» 



0.00510 



-0. 00041 
-0.00040 



0.0041 ±0.0011 
725 ± 22 



0.02779 ± 0.00020 [A08] 
0.253 ± 0.012 [A08] 

(fixed) [A08 & B08] 
(fixed) [A08 & B08] 
11.76 ±0.51 [B08] 
-1.48 ±0.93 [B08] 



678 ± 17 [A08] 



R* 
days 



'Transit epoch To 
'Orbital period P 
"Cl 

*c 2 



2454237.53556l^;^i 
1.7429935 ± 0.0000010 
0.911 ±0.016 
-0.094 ± 0.078 



2454237.53562 ± 0.00014 [A08] 
1.7429964 ± 0.0000017 [A08] 
0.88 ± 0.07 [A08] 
0.29 ± 0.07 [A08] 



HJD 

days 



Deduced parameters 



"Vsinl 


10.87 ±0.32 


11.85 ±0.50 [B08] 


km s _1 


P 


-4-0l 6 s i 


-7.2 ± 4.5 [B08] 


degrees 


RV K 


603 ± 18 


563 ± 14 [A08] 


rns -1 


btransit 


221 + 017 

u -" ± — 0.019 


0.253 ± 0.012 [A08] 


R* 


boccultation 


226+ 0018 


0.253 ± 0.012 [A08] 


R* 


Tpccultation 


2454238.40380lg;^o68 


2454238.40712 ± 0.00014 [A08] 


HJD 


Orbital semi-major axis a 


n n97QS+ 00076 

v.uz, l fO_ . 00080 


0.0281 ± 0.0009 [A08] 


AU 


Orbital inclination i 


88.081°:!! 


87.84 ± 0.1 [A08] 


degrees 


Orbital eccentricity e 


ni43 +0 0077 

U.UL43_ 0076 


[fixed, A08 & B08] 




Argument of pcriastron to 


10211, 7 




degrees 


Stellar mass M* 


0.96 ±0.08 


0.97 ± 0.06 [A08] 


M Q 


Stellar radius R* 


Qn6 +0 026 
u.yuo_ 027 


0.902 ± 0.018 [A08] 


Rq 


Stellar density p* 


1 000+0.035 
l.ZOO_o 033 


1.327 ± 0.006 [A08] 


p@ 


Wl 


346+ 014 

U.O^D_ oig 


0.41 ± 0.03 [A08] 






0.220 ± 0.032 


0.06 ± 0.03 [A08] 




Planet mass M p 


3.47 ±0.22 


3.31 ± 0.16 [A08] 


Mj 


Planet radius R p 


1 466+ 042 


1.465 ± 0.029 [A08] 


Rj 


Planet density p v 


1 inc+0.060 
1.1UO +0 056 


1.05 ± 0.08 [A08] 


Pj 



Table 1. CoRoT-2 system parameters and \-a error limits derived from our MCMC analysis. *A Bayesian penalty 
was used for the parameters preceded by an asterisk. 



measurement, extending thus considerably the time base- 
line. 

We thus confirm the low density of the massive planet 
CoRoT-2b. Tight constraints on the age of the sys- 
tem could help understand this peculiarity. As shown in 
Sec. 3.7, stellar evolution modelling does not much con- 
strain the age of CoRoT-2 (r* = 2.7l^ Gyr). Still, we 
have for this system three different age indicators. First, 
the presence of the Li I absorption line (B08) suggests a 



star still close to the ZAMS. We have obtained a new 
high-resolution high-SN spectrum of the star with the 
UVES spectrograph on the VLT (program 080.C-0661D, 
PI F. Bouchy). The Li I line at 6707 A is clearly de- 
tected in this spectrum (see Fig. 6). The Li I abundance 
measured from this line and MARCS atmospheric models 
(Gustafsson et al. 2008) is log n(Li) = ±2.8. Basing on 
Sestito & Randich (2005), this abundance suggests an age 
between 30 and 316 Myr. Secondly, the strong emission 
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line core in the Ca II H and K lines observed in the se- 
ries of HARPS spectra (B08) also indicates a young age. 
Using only the 14 HARPS spectra presenting a SN larger 
than 2 in the spectral region of the Ca II lines and the 
value 0.854 for the B — V color of the star from the Exo- 
Dat database (Deleuil et al. 2009), we deduce a value of 
-4.471 ± 0.0629 for the \og(R' HK ) parameter. Following 
Wright et al. (2004), this activity leads to an age of 307 
± 150 Myr. Finally, the small rotation period of the star 
measured from the CoRoT light curve (P ro t = 4.54 days, 
L09) and the B — V color of the star inserted into the 
relationship presented in Barnes (2007) lead to an age of 
76 ± 7 Myr for the system. As chromospheric activity is a 
magnetic phenomenon driven by rotation, it is clear that 
the \og(R' HK ) is not an age indicator independent of the 
rotation period (e.g. Vaughan et al. 1981, Noyes et al. 
1984). Still, we can conclude safely from both the large 
Li I abundance and rotational velocity of the star that 
CoRoT-2 is a very young star, of a few hundreds of Myr 
at most. In this context, a part of the 'radius excess' of 
the planet is explained. For instance, Fortney et al. (2007) 
models of irradiated planets predict a radius of ~ 1.3 Rj 
for a planet of 4 Mj orbiting at 0.02 AU around a 300 
Myr old solar-type star. 

Our occultation photometry imposes a strong con- 
straint on the parameter e cos ui and reveals that it is sig- 
nificantly smaller than zero. We can thus conclude that 
the orbit of CoRoT-2 is slightly eccentric. Unfortunately, 
e sin to is much less constrained by our data, so the ac- 
tual values of e and u> are poorly known, as shown by 
their marginalized posterior distribution function (PDF, 
see Fig. 7). The PDF of e itself is strongly not Gaussian: 
its 68.7% and 99.9% probability intervals are respectively 
0.007 < e < 0.022 and 0.001 < e < 0.037. To test the influ- 
ence of the Bayesian penalty on To and P on the resulting 
PDF of e cos w, we performed a new analysis without using 
these penalties. We obtained P = 1.7429926 ± 0.0000015 

and ecosw = -0.00258±o:ow)67- Thus the obtamed period 
does not disagree significantly (~ 1.7 sigma) with the one 
obtained by A08 from the CoRoT photometry, while the 
offset of the occultation remains significative (3.7 sigma, 
vs 4.6 sigma using the Bayesian penalties on P and To). 
This offset could be due to a slight eccentricity of the 
orbit, but also to a dynamical interaction with another 
object in the system (see, e.g., Schneider 2004, Holman & 
Murray 2005, Agol et al. 2005). Still, both the TTV anal- 
ysis presented in Alonso et al. (2009a) and the agreement 
between our deduced period and the one obtained by A08 
argue against this last hypothesis. We conclude thus to a 
slight eccentricity of the planetary orbit. 

Our result for the spin-orbit projected angle (3 agrees 
well with the one reported in B08, confirming a value close 
to zero for this parameter. To assess the influence of our 
Bayesian penalty on the parameter V sin /, we performed 
another MCMC integration without this penalty. We ob- 
tained Tasini = 11.6 ± 0.5 and — —5 ± 7 degrees, i.e. 
in agreement with the ones obtained with the Bayesian 
penalty, while becing slightly less precise. 




6708 6710 
Wavelength (A) 



Fig. 6. Portion of the CoRoT-2 UVES spectrum showing 
the Li I line and the weak contaminating Fc I line. 




0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 100 120 140 160 180 200 220 240 260 
e w (degree) 

Fig. 7. Marginalized PDF obtained for (from top left to 
bottom right) ecosw, esinw, e, and u>. Notice how the 
PDF for e and lo are non-Gaussian. 



Interestingly, the Spitzer flux estimator online tooJH 
indicates that the differences in magnitude that we mea- 
sured by deconvolution photometry between the nearby 
star and CoRoT-2 at 4.5 /zm (+1.7) and 8 ^m (+1.4) are 
consistent with a late-K or early-M type dwarf star located 
at the same distance (~ 200 pc) as CoRoT-2. As noticed in 
A08, this is also the case for the optical magnitude differ- 
ences from the Exo-dat database and the 2MASS near-IR 
magnitudes. In case of gravitational bounding, the angular 
distance between both stars would correspond to a phys- 
ical separation of ~ 800 AU. It is thus desirable to assess 
this possible gravitational bounding by independent mea- 
surements (proper motion, radial velocity). In case of con- 
firmation, CoRoT-2b would follow the tendency for mas- 
sive planets to be found preferentially in multiple stellar 
systems (Eggenberger et al. 2004). 
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4.2. Investigating the large radius of CoRoT-2b with 
coupled tidal-orbital evolution modelling 

CoRoT-2b is just one of many transiting planets with a 
radius larger than can be accomodated by standard ther- 
mal evolution models. Given the relatively young age of 
the planet, compared to other known transiting planets, 
it is worthwhile to investigate the planet's radius evolu- 
tion in some detail, as giant planets are expected to have 
larger radii at young ages. We use the coupled giant planet 
tidal and thermal evolution model of Miller et al. (2009) 
to calculate the planet's evolution and contraction. As in 
Miller et al., the planet's structure is assumed to have 
three components: a 50% rock 50% ice core, a fully convec- 
tive hydrogen-helium envelope with the equation of state 
of Saumon et al. (1995) and a non-grey atmosphere model 
described by Fortney et al. (2007) . The tidal orbital evolu- 
tion is described by Jackson et al. (2008a, 2009). This tidal 
evolution model assumes that the planet quickly reaches 
a spin-orbit sychronous state, that the only important 
source of tidal heating is due to orbital circularization, 
and the model is second order in eccentricity. 

In order to determine if tidal heating can explain the 
large radius of CoRoT-2b, for a variety of tidal quality 
factors Q' p , Q' s , and core masses, a grid over initial semi- 
major axis and eccentricity are evolved forward in time. 
We search for instances in each of these evolution histo- 
ries for which the semi-major axis, eccentricity, and radius 
are simultaneously within their error ranges. We choose 
to limit the age between 20 Myrs and 400 Myrs. We find 
that in cases when the Q' p value is too high (Q p of 10 6 or 
10 6 5 ) there is not sufficient dissipation inside the planet 
to achieve the observed radius. However, for the cases of 
Q' s = 10 5 -10 6 and Q' p < 10 5 5 all of the observed parame- 
ters can be explained as a transient event. Evolution histo- 
ries that closely agreed with the observed parameters are 
shown in Fig. 8 . The radius at optical wavelengths that 
the planet would be observed to have during the transit is 
shown in the upper left. The semi-major axis of the orbit 
is shown in the upper right. The ratio of input tidal power 
to net radiated power is shown in the lower left. The eccen- 
tricity is shown in the lower right. (See Miller et al. 2009 
for further details.) In each panel, the runs correspond to 
models with no core (black), 10 M® core (red), 30 M s 
core (blue). A "standard" run without tidal effects with 
a 10 M® core is in dotted cyan. The model without tidal 
heating clearly cannot explain the planet's larger radius, 
even given the young system age. This analysis suggests 
that if the Q p value is 10 5 5 or smaller, then it is possible 
to explain this large radius as a transient event at the last 
stage of orbital circularization. Under this scenario, the 
planet is spiralling inwards at high speed to its final tidal 
disruption, and the fast rotation of the star would not 
only be due to its young age but also to the high rate of 
angular momentum transfer from the planet's orbit. With 
such values for Q' s and Q p , the future lifetime of CoRoT- 
2b is 20 Myr at most, which is a short duration on an 
astronomical timescale but is still much larger than the 



Transit Radius Semimajor Axis 




time [Myrs] time [Myrs] 

Fig. 8. Possible tidal evolution histories for CoRoT-2. In 
these cases: Q' p — 10 5 5 and Q' s = 10 5 . For these curves we 
assume that the planet has no core (black), 10 M® core 
(red) and 30 M® core (blue). The cyan run assumes that 
the planet has a 10 M® core with no tidal evolution. See 
text for discussion 

remaining lifetime of the planet WASP- 18b (Hellier et al. 
2009) under similar assumptions. 

In some planetary systems, an outer companion might 
continously drive the eccentricity of the inner planet off- 
seting circularization by tides such that the eccentricy is 
found in a semi-equilibrium state, described by Mardling 
(2007). Let us assume this scenario is occuring and the 
planet's net radiated luminosity, L p , at the surface is bal- 
anced by tidal heating inside, Ft- Using Table 1 from 
Miller et al. (2009) 

L p = 7 x 10 28 ergs/sec 




and assuming that the observed eccentricity of 0.0142 is 
close to its equilibrium value, then this would imply that 
Q' p ~ 10 4 . This is lower than the oft-quoted value for 
Jupiter between 10 5 and 10 6 (Goldreich & Soter, 1966). 
In summary, we find that a young age alone cannot ex- 
plain the large radius of CoRoT-2b, but that plausible 
tidal heating evolutionary histories, with Q' ~ 10 4 — 10 5 5 , 
can explain it. 

4.3. Atmospheric properties of the young planet 
CoRoT-2b 

The thermal emission of the planet is detected in 
both IRAC channels, as can be seen in Fig. 4 and 9. 
Unfortunately, the precision on the occultation depth at 
8 /im is rather poor and thus brings a weak constraint 
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on the planetary SED. CoRoT-2 is indeed a faint tar- 
get for Spitzer at 8 /im, the theoretical error (photon, 
read-out and background noise) per 10.4 s exposure being 
~ 0.84%, while it is ~ 0.29% at 4.5 fim. The standard 
deviation of the residuals of our best-fitting solution are 
close to these values: 1.09% (8 fxm) and 0.32 % (4.5 /jm), 
i.e., respectively, 1.3 and 1.1 times larger than the the- 
oretical noise budget. Assuming that the observed noise 
is purely white and taking into account the error on the 
flux normalization, we would expect an error of ~ 0.06% 
on the occultation depth at 8/im, while our MCMC anal- 
ysis, which takes into account the low-frequency noises, 
leads to an error ~ 1.8 times larger. We can thus conclude 
that the large level of correlated noise (due to the ramp, 
the low-frequency stellar and background variability, the 
blend with the nearby fainter star, and other unknown ef- 
fects) has a significant effect on the final precision. This is 
also the case at 4.5 /zm: for purely white noise, we would 
have expected a precision of 0.016% on the occultation 
depth, while our actual error is ~ 2.5 times larger. 

One of the many interesting questions that has arisen 
since the direct detection of hot- Jupiter atmospheres be- 
gan (Charbonneau et al. 2005; Deming et al. 2005) is 
which, if any, of these planets has temperature inversions 
in their atmospheres (see Sec. 1). So far, of the seven plan- 
ets with published Spitzer measurements in each of the 
four IRAC channels, all but HD 189733b and TrES-3b 
have been reported to exhibit temperature inversions (see 
references in Sec. 1). Figure 9 compares the IRAC colorsjf] 
for these seven hot Jupiters. The lack of a clear pattern in 
Fig. 10 illustrates the difficulty in using Spitzer photom- 
etry alone to identify an atmospheric inversion and high- 
lights the model dependencies of the inversions inferred 
so far from these data. The scatter in this diagram also 
demonstrates that an observationally based classification 
scheme cannot yet be defined. 

For CoRoT-2b, optical measurements are also available 
(Alonso et al. 2009b, Snellen et al. 2009b). These latters 
and our two IRAC measurements are compared in Fig. 11 
and Table 2 to three different models: 

— Model 1 (mi) assuming an efficient heat distribution 
to the night-side of the planet. 

— Model 2 (7712) assuming no heat distribution to the 
night-side and no temperature inversion. 

— Model 3 (7/2.3) assuming no heat distribution to the 
night-side and a deep TiO/VO-induced temperature 
inversion. 

Table 2 shows also the x 2 obtained for each model. 
Given the \ 2 of two models m a and nib, we can compute 
their likelihood ratio: 

(X 2 (>na)-X 2 (»n b )) 

LR(m a /m b ) = e * (9) 

Comparing the model mi to the two others, we obtain 
LR{m 1 /m 2 ) = 3x lO" 11 and LR( mi /m 3 ) = 2 x 10~ 9 . An 

6 Here color is calculated by taking the ratio of the planet 
fluxes in the IRAC channels and have not been scaled by the 
flux ratios for Vega. 
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Fig. 9. Marginalized PDF for the IRAC occultation depth 
at 4.5 itm (top) and 8 /im (bottom). 
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Fig. 10. IRAC color-color diagram for six hot Jupiters. 
Symbol sizes are scaled by the level of incident starlight 
received by the planet. The Fortney et al. (2008) classifi- 
cation is also indicated. The location of CoroT-2 falls in 
the shade region, indicating the l-cr uncertaities for the 
[4.5]-[8.0] color. The solid arrow indicates the trend for 
blackbody-emitting planets. See text for references. 



efficient heat distribution to the night-side of the planet 
is thus strongly disfavored by the data. Comparing now 
the models 7712 and 7773, we obtain LR(m,2/'m 3 ) — 74. The 
data are thus in better agreement with the absence of 
a strong temperature inversion, but the obtained likeli- 
hood ratio is not large enough to firmy conclude, and a 
greater wavelength coverage and a better precision in the 
observations are required to remove the current ambiguity. 
Snellen et al. (2009b) concluded too that the C0R0T opti- 
cal measurements favor an absence of temperature inver- 
sion. Nevertheless, we notice that their best-fitting model 
significantly under predicts the flux measured at 4.5 /xm, 
this latter being in much better agreement with their mod- 
els assuming an inversion. 
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CoRoT white 


CoRoT red 


IRAC 4.5 pm 


IRAC 8 pm 


x 2 


Measured 


0.0060 ± 0.0020* 


0.0102 ± 0.0020" 


0.510 ±0.041 


0.41 ±0.11 




Model 1 


0.0021 


0.0022 


0.270 


0.418 


55.1 


Model 2 


0.0062 


0.0075 


0.449 


0.587 


6.6 


Model 3 


0.0068 


0.0053 


0.528 


0.740 


15.2 



Table 2. Comparison of the measured planet-to-star flux ratio and the values predicted by our three models. See text 
for details. *from Alonso et al. (2009b). **from Snellen et al. (2009b). 
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Fig. 11. Synthetic planet-star flux density ratios from 
three hot-Jupiter atmosphere models (from Barman et al. 
2005). The top two curves are from models with incident 
stellar flux constrained to the dayside, while the lower 
curve corresponds to a model with unifom day-to-night 
redistribution of stellar flux. The dotted line shows the 
flux-ratio for a planet with a TiO/ VO-induced tempera- 
ture inversion. Solid points, with l-cr error-bars, are the 
Spitzer-IRAC and CoRoT photometry (see inset). Open 
symbols are the band- integrated model points. 

5. Conclusion 

Using Spitzer and its IRAC camera, we measured an oc- 
cultation of the young and massive planet CoRoT-2b at 
4.5 /jm and 8 pm. In addition, we observed a partial tran- 
sit of the planet with the Very Large Telescope and its 
FORS2 camera. 

Our global analysis of CoRoT, Spitzer and ground- 
based (FORS2 photometry ± RVs) data confirms the 
low density of the planet (p p = 1.105 ± 0.060 pj with 
M p = 3.47 ± 0.22 Mj and R p = 1.466 ± 0.044 Rj). 
Constraining the system to be of at most a few hundreds 
of Myrs old and the present orbit to be slightly eccentric 
(ecosw = -0.00291 ± 0.00063) and using coupled tidal- 
orbital evolution modelling, we find a self consistent ther- 
mal & tidal evolution history that may explain the radius 



through a transient tidal circularization & corresponding 
tidal heating inside the interior of the planet. Under this 
scenario, the planet will be tidally disrupted within 20 
Myr at most. 

The occultation depths that we measured at 4.5 //m 
and 8 ^m are, respectively, 0.510±0.042 % and 0.41 ±0.11 
%. In addition to the optical measurements reported by 
Alonso et al. (2009b) and Snellen et al. (2009b), these 
values favor a poor heat distribution to the night side of 
the planet and the absence of thermal inversion, but mea- 
surements at other wavelengths are needed to confirm this 
latter point. 
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